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Abstract 

Asymptotic formulae for Green's kernels G£(x,y) of various boundary value prob- 
lems for the Laplace operator are obtained in regularly perturbed domains and 
certain domains with small singular perturbations of the boundary, as e ^ 0. The 
main new feature of these asymptotic formulae is their uniformity with respect to 
the independent variables x and y. 



Introduction. In 1907, Hadamard was awarded the Prix Vaillant by tlie 
Academie des Sciences de Paris for his work [1]. Among much else, this mem- 
oir includes investigation of Green's kernels of some boundary value problems, 
and, in particular, their dependence on a small regular perturbation of the do- 
main. Hadamard's asymptotic formulae played a significant role at the dawn 
of functional analysis (see [2]). The list of subsequent applications includes 
extremal problems in the complex function theory [3] , [4] , a biharmonic max- 
imum principle for hyperbolic surfaces [5], shape sensitivity and optimisation 
analysis [7], free boundary problems [8], Brownian motion on hypersurfaces 
[6], theory of reproducing kernels [9], [10], [11]. Analogues of Hadamard's for- 
mulae were obtained for general elliptic boundary value problems [13] and for 
the heat equation [14]. 



The issue of asymptotic forniTilae for Green's kernels in singularly perturbed 
domains attracted less attention. In this respect, we mention the papers [15], 
[16], where certain estimates for Green's functions in domains with small holes 
were obtained. 

An important drawback of the estimates just mentioned, which is also inher- 
ent in the classical Hadamard formulae, is the non-uniformity with respect 
to the independent variables. The present article settles the question of uni- 
formity left open in earlier work, both for regularly and singularly perturbed 
domains. First, we sharpen one of Hadamard's classical formulae in the case 
of a regular perturbation. Next, we consider several types of singularly per- 
turbed domains including domains with a finite number of holes, thin rods and 
a truncated cone. For these geometries and different types of boundary condi- 
tions, we derive uniform asymptotic representations of Green's kernels for the 
Laplacian. The asymptotic formulae presented here for singularly perturbed 
problems involve Green's functions and other auxiliary solutions of certain 
model boundary value problems independent of e. 

The techniques used in the proofs of theorems of the main text of the paper 
are based on the method of compound asymptotic expansions discussed in 
detail in the two- volume monograph [17]. The uniform asymptotic formulae 
for Green's kernels can be efficiently used in numerical algorithms for solving 
boundary value problems in domains with singularly perturbed boundaries. 

1. Uniform Hadamard's type formula. Let be a planar domain with 
compact closure Q and smooth boundary Also, let another domain Jl(£), 
depending on a small positive parameter e, lie inside Vt. By 5z we denote a 
smooth positive function defined on dVt^ and assume that e5z is the distance 
between a point z e dVL and dVte- One of the results in Hadamard's paper 
[1] is the following formula, which relates Green's functions G and for the 
Dirichlet boundary value problem for the Laplacian in Q and Q.^'. 



This asymptotic relation holds, in particular, when either x or y is placed at 
a positive distance from dVL, independent of e. However, one can see from the 
simplest example of two concentric disks and Jl^ that (1) may fail when 
x and y approach the same point at dO.^. In the next theorem, we improve 
Hadamard's formula by obtaining a uniform asymptotic representation for 



Theorem 1. Let x and y he points ofQ^ situated in a sufficiently thin neigh- 



Ge(x,y)-G(x,y)+£ 





G,-G. 
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bourhood of dD,. The right-hand side in (1) can be replaced by 
1 , - z(y)p + (p^ - eS,^^) +Py- 

An ^ |^(x)-;,(y)|2 + (p^+p^)2 ^ ) 



^27r |^(x) - z{y)\' + (p, + p,)^ ^ ^|^(x) - z{y)\' + (p. + p,)^^' 

where z(x) t/ie pomi o/9Q nearest to x, p^: = |x — z(x)|, anc? ^/le notations 
z{y), py have the same meaning for the point y. 

We note that the principal term in (2) plays the role of a boundary layer, and 
when |z(x) — z(y)| < Const e and Px + Py < Const £, it has the order 
whereas the remainder term has the order 0{e). 

2. Dirichlet problem for a domain with a small inclusion. Let Q and 

uj be bounded domains in R". We assume that Q and u contain the origin O 
and introduce the domain a;£ = {x : e^^x G ur}. Without loss of generality, 
it is assumed that the minimum distance between the origin and the points 
of dQ as well as the maximum distance between the origin and the points 
of duj are equal to 1. Let be Green's function of the Dirichlet problem 
for the Laplace operator in fi^ = \ cJ^. We use the notation l^*""^! for the 
(n — l)-dimensional measure of the unit sphere. 

Theorem 2. Let n > 2. By G and <3 we denote Green's functions of the 
Dirichelt problems in VL and R"' \ W, respectively. 

LetH be the regular part ofG, that is H{x,y) = ((n-2)|5"-i|)-i|x-yp-"- 
^(x, y), and let P stand for the harmonic capacitary potential ofuJ. Then 

G,{x, y) = ^(x, y) + £2-0(£-ix, £" V) - {{n - 2)\S^-'\)-'\x - y|2- 
+H{0, y)P(£-^x) + i/(x, 0)P(£- V) - H{0, 0)P(£-ix)P(£- V) 

- £"-2 cap uJ H{^, Q)H{0, y) + o( . ' ) (3) 

^(mm{|x|, |y|})" 



uniformly with respect to x and y in^.^. (Note that the remainder term in (9) 
is 0{e) on du^ and 0{e'^~^) on dQ. ) 

The next theorem contains a result of the same nature for n = 2. As before, 
G and (S are Green's functions for ^2 and R^ \ LJ, respectively, whereas H is 
the regular part of G. 
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Theorem 3. Let 

C(t7)= lim 05(^,77) and Coo = Hm {((r?) - (27r)-Mog |r;|}. 

Then the asymptotic representation, uniform with respect to x,-y E fls, holds 
G,(x, y) = G(x, y) + 0(£-^x, e'V) + log(£-i|x - y|) 

((27r)-^ logg + C(f ) - Coo + Hj^, 0)) ((27r)-i loge + C(f ) - Coo + HjO, y)) 
^ (27r)-Mog£ + i/(0,0)-Coo 

-C(£-'x)-C(£-V) + Coo + 0(£). (4) 



The next assertion is a direct consequence of (4). It shows that asymptotic 
representation of ^^(x, y) is simplified if x and y are subject to additional 
constraints. We use the same notations as in Theorem 3. 

Corollary. The following assertions hold: 

(a) Let X and y be points of fl^ C R"" such that min{|x|, |y|} > 2e. Then for 
n = 2, 

r( \ r( \ G'(x,0)G(0,y) e 
G£(x,y) - G'(x,y) - — — — , ^ = ^'l 



(27r)-i log£ + H{Q, 0) - Coo 'min{|x|, |y|}' 
and for n > 2, 

-n-l 



G,{^, y) - ^(x, y) + £"-2cap U ^(x, 0)^(0, y) = 0( ^ 



(b) If max{|x|, |y|} < 1/2, then for n = 2, 



(|x||y|)" imin{|x|, |y|}' 



and for n > 2, 



G,(x,y)-£2-"0(-,^)+i/(O,O)(P(£-^x)-l)(P(£-V)-l) = O(max{|x|,|y|}). 

e e 



One can say that these formulae are closer, in their spirit, to Hadamard's 
original formula (1) than Theorems 2 and 3 themselves. Analogous simplified 
formulae can be deduced easily from all other theorems of the present paper, 
but they are not included here. 
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3. Dirichlet-Neumann problems in a planar domain with a small 
hole. The set Q,, is assumed to be the same as in Theorem 3, with the ad- 
ditional constraint that duj is smooth. First, let denote the kernel of the 
inverse operator of the mixed boundary value problem in for the operator 
—A, with the Dirichlet data on dQ and the Neumann data on dcoe- The no- 
tations G, ^5 and H have the same meaning as in Theorem 3, and 91 is the 
Neumann function in \ uJ. Let 2) be a vector function, harmonic in \ uJ, 
vanishing at infinity and such that d'J^/du = u on du>. This vector function 
appears in the asymptotic representation 

9t(e»7)~(27r)-Mog|£|-i + (5:)(r7)-77)-V((27r)-Mog|£|-i), as |£H oo. 



Theorem 4. The kernel (je(x, y) of the inverse operator of the mixed boundary 
value problem in fi^ for the operator —A, with the Dirichlet data on dO, and 
the Neumann data on du^, has the following asymptotic representation, which 
is uniform with respect to x, y e fig, 

^.(x, y) = G(x, y) + OT(£-^x, e-^y) + (27r)-^ log(£-^|x - y|) 

+£S)(£-^x) • V,i/(0, y) + V) • V,i/(x, 0) 

- £^5?(£-V) • ((V, V,)//(0, 0))2)(£-ix) + 0{e^). (5) 



Next, consider the mixed boundary value problem in Vt^ for the Laplace oper- 
ator, with the Neumann data on the smooth boundary dVL and the Dirichlet 
data on duj^, where a;^ = {x : £~^x e a; C R^} with uj being an arbitray 
bounded domain. Let A'"(x, y) be the Neumann function in Q, i.e. 

AAr(x,y)+(^(x-y) = 0, x, y e O, 

d 



■(7V(x,y) + (27r)-Mog|x|) =0, ^edn,yen, 



and 

/d 
N{x, y)— log|x|ds^ = 0, 

on ^ 

with the last condition implying the symmetry of N{x.,y). The regular part 
of the Neumann function is defined by 

i?(x, y) = (27r)-^ log |x - y]"^ - 7V(x, y). 

Note that 

R{Q,y) = -(27r)-2y' log |x| — log |x|ds^. 

an 
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Let D be a vector function, harmonic in R^\a7, bounded at infinity and subject 
to the Dirichlet condition D(^) = ^, ^ G duo. This vector function appears 
in the asymptotic representation of Green's function in \ UJ 

0(^,77)~0(oo,77) + (D(r7)-77)-V((27r)-Mog|£|-i), as |£H oo. 

Theorem 5. The function Gei'x., y) has the asymptotic behaviour, uniform 
with respect to x, y e Q^: 

G,(x, y) = 0(5-ix, e-'y) + iV(x, y) - {2n)-' log |x - y\-' + i?(0, 0) 

+£D(£-V) • V^i?(x, 0) + £D(£-^x) • V^i?(0, y) 

- s'B{e-'y) ■ (((V. ® V,)i?)(0, 0))d(£-1x) + 0(5^). (6) 

4. The Dirichlet-Neumann problem in a thin rod. Let C be the infinite 
cylinder {(x',a;„) : x' G cj, a;„ G R}, where is a bounded domain in 
R"~i with smooth boundary. Also let denote Lipschitz subdomains of 
C separated from ±oo by surfaces 7^, respectively. We introduce a positive 
number a, the vector a = (0',a), where O' is the origin of R""-*^, and the 
small parameter £ > 0,. 

By Ge we denote the fundamental solution for —A in the thin domain = 
{x : £~^(x — a) e C+, £'~^(x + a) G C"} subject to zero Neumann condition 

on the cylindrical part of C and zero Dirichlet condition on the remaining part 
of dCg. Similarly, and G^ stand for the fundamental solutions for —A in 
the domains C^, and satisfy zero Dirichlet condition on 7^, zero Neumann 
condition on dG^ \ 7^, and are bounded as Xn ^ Too. Let be harmonic 
functions in subject to the same boundary conditions as and asymptot- 
ically equivalent to =F|^|~''"2;n + Too, where arc constants. We 
note that C^(y) = li^ix^oo y). Finally, by Goo(x, y) we denote Green's 
function of the Neumann problem in G such that 

Goo(x,y) = -(2|a;|)~^|x„ +0(exp(-Q;|x„ -y„|)) as ^ 00, 

where a is a positive constant, and |ci;| is the (n — l)-dimensional measure of 
cu. 

Theorem 6. The following asymptotic formula for Gs{'x,y), uniform with 
respect to x, y G fie, holds 

G,(x,y) = £2-|G+(£-Hx-a),£-Hy-a))+G-(£-i(x+a),£-i(y+a))-G-(£-ix,£-V) 

^ i^l - l(C - C+) + C+(^) - C-(^))(fe - hiC^ - C+) + C+(^) - c-(^)) 

2|u;|-ia + e(Ci) + C00) 
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+\ {{e\u\r'2a + C- + C+ - 2 (C^(£-'(x T a)) + (^(^-^(y t a)))) 

f 0(exp(-/?/£))}, (7) 



where (3 is a positive constant independent of e. 

5. The Dirichlet problem in a truncated cone. Let K be an infinite cone 
{x : |x| > 0, |x|^^x G tu}, where a; is a subdomain of the unit sphere S'^"^ 
such that S'"~^ \ u; has a positive (n — 1) -dimensional harmonic capacity. Also 
let Kq and K^o denote subdomains of K separated from the vertex of K and 
from oo by surfaces 7 and F, respectively. 

By Ge and Gcone we denote Green's function of the Dirichlet problem for —A 
in the domains = {x & Kq : e~^x G K^o} and K, respectively. Similarly, 
Go and Goc stand for Green's functions of the Dirichlet problem for —A in the 
domains Kq an K^o satisfying zero Dirichlet boundary conditions on dKo\{0} 
and dKoo, with the asymptotic representations 

Zo(x) = |x|2-"-^*(|x|-^x)(l + 0(1)) as |x| ^ 0, 

and 

Zoo(x) = |x|^*(|x|~^x)(l + 0(1)) as |x| 00, 

where A is a positive number such that A(A + n — 2) is the first eigenvalue of 
the Dirichlet spectral problem in u for the Beltrami operator on 5"^"^, and \E' 
is the corresponding eigenfunction. Similarly, we introduce A2 > such that 
A2(A2 + n — 2) is the second eigenvalue of the same spectral problem. 

Theorem 7. Let £ = £~^x and rj — e~^y. Green's function ^^(x, y) has the 
asymptotic behaviour 

G,(x, y) = Go(x, y) + e^'^G^it rj) - G^^nei^, y) 

+ {\ri\'^{^) - ^oo(r7))(|x|2-"-^v]>(A) _ Zo(x))} + 0(£--(2^'^^)).(8) 
This representation is uniform with respect to x and y in the truncated cone 

6. Dirichlet problem in a domain containing several inclusions. It 

is straightforward to generalise the results of sections dealing with a domain 
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containing a small inclusion/void to the case of a body containing a finite 
number of inclusions. As an example, we formulate a generalisation of Theorem 
2. Let C R",n > 2, be a bounded domain, and let 0^^\0^^\ O^^) be 
interior points in Q. Small sets LO^J\j = 1, . . . , A^, are defined by uj^J^ = {x : 
e~^{yi — O^^) e ijj^^ C R"}, where uj^\ j — 1, . . . ,N, are bounded domains 
in R", and they contain the origin O. Similar to Section 2, it is assumed that 
the minimum distance between 0'^^\j = 1,...,N, and the points of dVl is 
equal to 1. Also, it is supposed that the distance between O and the points of 
duj^^\j = 1, . . . , A", does not exceed 1. By we denote Green's function for 
the Laplacian in the domain D,^ — Q,\ [JjUJ^^\ 

Theorem 8. Let G and ^^'-i^ stand for Greenes functions of the Dirichlet prob- 
lems in VL and R" \ respectively. Also let H he the regular part of G, and 
P^^^ denote the harmonic capacitary potentials of the sets cJ^\ Then 

^ , ^ X - 0(-^') V - O^^') A" 



j=i - - ,n-2)|5^-i||x-y 
+ E {H{0^^\y)P^'\^^ ) + H{^, 0(^))p(^\l^^ ) (9) 



n-2 



3=1 



H(0(^\ QO))pO)( '' capu;(^) HU, O^^AH(0^^\y)] 

£ £ ' 



( (min{|x-OW|,|y-O0')|})"-2) 

uniformly with respect to x and y infl^. (Note that the remainder term in (9) 
is 0{£) on 9(Uja;p)) and 0(£"-^) on dfl. ) 
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